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Abstract 

We consider the single particle relaxation rate of 2D electrons subject to a 

random magnetic field. The density of states (DOS) oscillations in a uniform 

magnetic field (in addition to a random one) are studied. It is shown that the 

infrared singularities appearing in perturbation theory for the single particle 

relaxation rate are spurious and can be regularized by a proper choice of gauge. 

We define a gauge invariant single particle relaxation time as a parameter 

describing the amplitude of DOS oscillations. Unlike the conventional case 

of random potential scattering the broadening of the Landau levels does not 

depend on magnetic field. 

PACS numbers:72.15.Lh, 71.25 He 
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The problem of a charged quantum particle moving in a 2D static random magnetic field 
has attracted considerable interest recently. Models of this type are supposed to describe 
effectively a number of important physical systems, such as a state with spin-charge sepa- 
ration in high-T c superconductors and a quantum Hall system near the filling factor 
v = 1/2 j|,|5||. In addition, the direct experimental realization of static random magnetic 
field has been reported recently |J. 

In physical realizations of such a system one is interested in a random magnetic field 
distribution characterized by some correlation length. This is unlike the usual case of a 
disordered system, where the electric potential is a random variable with finite range corre- 
lations. Here the corresponding quantity is the vector potential, which will be governed by a 
spatially long-ranged correlation function, even for short-ranged magnetic field correlations. 
Consequently, similar to the case of long-range potential fluctuations, small angle scatter- 
ing will play an important role, and might even cause the total scattering cross-section to 
diverge. On the other hand, the transport scattering cross-section, which is weighted with 
the momentum transfer squared, is expected to be well-behaved. 

In a previous paper J7| the transport properties of the model were studied. It was shown 
there that although the single particle relaxation rate 1/r is plagued by infrared divergence 
problems, the electrical conductivity is determined by the transport (momentum) relaxation 
time Tt r , which can be calculated perturbatively (see also |§) . As far as the critical behavior 
of the model is concerned, the problem could be mapped onto a nonlinear supersymmetric a- 
model of interacting matrices with unitary symmetry. The diffusion constant (or equivalently 
the dimensionless conductance) go was shown to be determined by r tr rather than r. This 
means that all the states are localized within the localization length £ oc exp (irgo) 2 . On 
shorter length scales the transport is governed by the usual diffusion law, in contrast to the 
suggestion made in || that the particle diffuses according to a logarithmic diffusive law. 
We believe this to be an artefact of an improper approximation employed in 0, as will be 
discussed in the end of the paper. 

In the present paper we consider the single particle properties of the model. There are 
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two related problems that have to be dealt with: the existence of infrared divergences in 
perturbation theory and the gauge dependence of the Green's function. The latter has been 
studied recently |TIJ in a path integral formulation. In order to circumvent the gauge problem 
the authors of JIIJ proposed to introduce a gauge invariant Green's function, obtained by 
multiplying the usual single particle Green's function by a phase factor. The phase was 
chosen as the one caused by integrating the vector potential along a straight line connecting 
the end points. The Green's function defined in this way was evaluated in the quasiclassical 
approximation. Obviously, this quantity depends on the choice of the integration path 
determining the phase. The physical meaning and possible applications of this Green's 



function are far from clear. The same object later was studied in ||, and the results of flO 
were reproduced. 

We define below the single particle relaxation time as the parameter of the broadening 
of the Landau levels in a weak, constant magnetic field (in the presence of the random one). 
This is an unambiguous experimentally meaningful parameter which is definitely gauge 
invariant. It can be determined by measuring the de Haas - van Alp hen or Shubnikov-de 
Haas oscillations. Note that the constant magnetic field appears naturally and inevitably 
for the direct experimental realization of the model ||. 

In a previous paper @ it was observed that the calculation of the single particle relaxation 
time is complicated by the infrared divergence of the corresponding diagrams. Meeting 
analogous problems in the slave-boson representation of the t — j model the authors of 
concluded that the two-point auxiliary particle Green's function should vanish identically. 
These authors concentrated on the consideration of the transport time. We will show here 
that the infrared divergence is fictitious and can be removed by a gauge transformation. 
The procedure we use for this purpose is somewhat analogous to that used in PHJI3]] for 
separation of the low-momentum transfer scattering processes in a random potential with 
long-ranged correlations. Although the resulting Green's function has no infrared singularity 
it still changes under gauge transformations. This should not pose a problem for the one 
point Green's function, which is gauge invariant and the imaginary part of which determines 



the density of states (DOS). 



We show that the use of the self-consistent Born approximation p?|-|T4|j for calculation 
of the self-energy violates the gauge invariance and leads to an ambiguous result for r. 
The ambiguity can be resolved, however, by using the formalism of path integrals |15j in 
coordinate space. We obtain a Gaussian, rather than Lorentzian, shape of the broadening of 
the Landau levels which explains the difficulties arising when one tries to treat the problem 
by means of the conventional diagram technique. 

Consider a charged spinless particle (with the mass m and charge e) in a 2D layer, in 
a random magnetic field H(r) perpendicular to the layer. The random magnetic field is 
supposed to be Gaussian distributed with a zero mean value and correlator: 

(H(r)H(r')} = (h 2 )5(r - r') (1) 

The vector potential correlator, in momentum space, can be taken in the form: 

(Aa(q)M-Q)) = <W^r = ( 2 ) 

Applying the standard perturbation theory with the correlation function (fj), one encounters 
divergences generated by the singular nature of (Q) in the limit q — > 0. A similar problem 
was considered in [llJHj for the case of the long-ranged Coulomb potential. Exploiting the 



ideas from [11,13|] we represent the correlator D(q) in the form: 

D(q)=D (q)+D{q) 

D {q) = d 8(q); d = f d 2 qD(q) (3) 

Jg<Q 

where Q is an arbitrary momentum cut-off, which serves to separate low- and high- mo- 
mentum contributions. It can be shown that the averaged single-particle Green's function 
G(E,p) = / E — 2^;(p — f^4-) 2 ) can be rewritten in the form : 



where G differs from G in that the correlator D(q) (rather than D(q)) is used for the 
averaging over the disorder. 



The crucial difference between (HI) and the analogous expression in |TT|,p!3| is that on 



the r.h.s. of (j3]) the averaging enters by the shift of the momentum p, whereas in pH|J13 
it entered as a shift of the chemical potential. In principle, the transformation (Q) affects 
the Green's function in a non-trivial way. However, if we consider the one-point Green's 
function / d 2 p G(E,p) |TjJ, shifting the integration variable p we obtain: 

Jd 2 pG(E,p)= Jd 2 pG(E,p) (5) 

This is nothing else than the manifestation of gauge invariance. The correlator D (q) corre- 
sponds to a spatially uniform vector potential, i. e. to zero magnetic field. Thus, one can use 
D(q) instead of D(q) for the calculation, since it does not affect any gauge invariant quantity. 
The advantage of this simple trick is that now all the infrared divergences are removed from 
the diagrams. This can be checked for any particular diagram. Demonstrate the cancelation 
of the infrared divergences for the sum of non-crossing diagrams, which represents SCBA. 
For simplicity we consider the case when the constant magnetic field is equal zero. We have 
for the imaginary part of the self-energy: 

d 2 q \ (2p + q) 2 fi/2T' (2p) 2 h/2r } 



h/2r(E,p)=v 2 f 

J a 



q <Q q 2 [(E — (p + q) 2 /2m) + (h/2r') (E - p 2 /2m) z + (ft/2rT J 
a r fq (2p+q) 2 h/2r' 
V ° J q >Q q 2 (E-(p+ q) 2 /2mf + {h/2r'f 1 ' 

where v 2 = e 2 (h 2 ) /4m 2 c 2 is the characteristic velocity of the problem, r = r(E,p), t' = 

r(E,p + q). It is easy to see that the first integral does not diverge at small q anymore. 

Since we are interested in the mass-shell value of r, we take E = p 2 /2m and suppose that 

Et ^> h. The integrals are dominated by a vicinity of the mass-shell, where one can take 

r' = t. A straightforward calculation gives : 

As was discussed above, due to the gauge symmetry, no observable quantity can depend 
on the cut-off parameter Q. However, repeating the same derivation for the case of non- 
zero constant magnetic field, one obtains the same ambivalent logarithm for the amplitude 



of oscillations of the density of states, which is definitely a gauge invariant quantity. The 
reason for obtaining such an ill-behaved single particle time is that we implicitly expected 
the Green's function to be of the usual, Lorentzian, form: 



Actually one should not expect this form of the Green's function. Assume, for a moment, 
such a form for a particular value of Q. Then any change of Q transforms the Green's 
function in a way described by equation (§]), and the Green's function acquires a Gaussian 
form rather than the Lorentzian one. Thus the assumption that the Green's function can 
be described by (§) contradicts the gauge invariance, and leads to the above mentioned ill 
behavior of the single particle time. 

The physical reason for the appearance of the non-gauge-invariant logarithm in (0) is 
the following. Consider the open trajectories for the two-point Green's function in a path 
integral representation. Gauge freedom means that we can add arbitrary far remote magnetic 
flux tubes to the system which will give rise to additional random phases. The choice of Q 
specifies the scale for the influence of the remote random fluxes on the Green's function. This 
ambiguity is removed if we restrict the consideration only to closed trajectories {i.e. one- 
point Green's function) on which the random phase is a well defined quantity determined 
by the random flux enclosed by the trajectory. However, in order to calculate the self- 
consistent single particle relaxation rate we are forced to consider open trajectories, which, 
having gauge ambiguous random phases, cause the appearance of the ambiguous logarithm 
in the corresponding relaxation time (0). 

To overcome these difficulties, arising in the diagrammatic approximation, let us consider 
a different approach to the problem. In the following we use the quasiclassical approximation 
i.e. assume Er ^> h and E/uj ^> h, where E is the Fermi energy, u = eH c /mc is the 
cyclotron frequency, H c is a uniform magnetic field. 

It is known |T5[ that the one-point Green's function can be obtained, in the quasiclassical 
limit, by calculating the path integral over the closed classical trajectories. In our case, after 




(8) 
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averaging over the random magnetic field, the action on such trajectories is 

2 (h 2 ) 



i-H c s D 



(9) 



Here s Q and s no denote the oriented and non-oriented (Amperean) areas of the trajectory 
respectively. The first term describes the phase acquired by the particle due to the flux of 
the constant magnetic field and the second appears due to averaging over the random one. 
For trajectories without self- intersections the areas can differ only in sign (s no = \s a \). The 
non-oriented area of a self-intersecting trajectory depends non-trivially on its topology, e.g. 
for the closed trajectory with a winding number k the non-oriented area is proportional to 
k 2 . 

We will consider the classical trajectories determined by the first term in ([]), and treat the 
second term as a perturbation. The criterion of the applicability of this approximation will be 
discussed later. Classical trajectories in a magnetic field are circles with a cyclotron radius 
R c = p/(mu). Calculating the action on the classical trajectory S(E) = (e/c)H c irR 2 c = 
2ttE/uj, we can write the Green's function as a sum over classical trajectories with different 



winding numbers [17]: 



G(E) 



—im 



1 + 2 ]T exp \ 2nki 



k=l 



E 1 
2 



Auk 2 - rffl!)? 

nuj 1 



(10) 



Comparing the two terms in the exponent in ( |10|) one can see that the perturbative treatment 
of the random field is valid, provided E 3> lo ^> itivq. To obtain explicitly the shape of the 
broadened Landau level we use the Poisson formula to resum (|i~0| ) and represent the DOS 
p(E) = — - ImG(E) as a sum over the energy quantum numbers N: 

1 



p(E) 



E - hu){N + - 
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(ii) 



where Ir = \Jch/eH c is the magnetic length. This formula evidently shows that the Landau 
levels acquire the Gaussian form. This fact explains the difficulties which one meets trying to 
use the conventional diagram technique, which implicitly assumes a Lorentzian broadening 
of the Landau levels. The Gaussian shape of the Landau levels appears here as a result 



of the quasiclassical approximation. For the case of the potential scattering one can carry 
out the analogous derivation, when the potential is smooth enough to use the quasiclassical 
approximation. This also leads to a Gaussian shape of the DOS in agreement with [JT3| . 
The characteristic width of the levels is : 

n n 2 \ 1/2 , , 

-Emvl) (12) 



2t(E) Vtt 

Note that this expression surprisingly coincides with one obtained by SCBA, except for the 
ambivalent logarithm. The cyclotron frequency does not enter the expression for r, thus one 
can consider r as a meaningful parameter even in the limit of zero magnetic field H c . 

Equation (|10D describes DOS as a function of energy. For the case of weak uniform 
magnetic field (lot < 1), omitting all the harmonics except the first one we obtain for the 
oscillating DOS p osc (E): 

^ (B) = __^ exp {__J_} cos ( 2 ^) (13) 

One can see that the amplitude of the DOS oscillations is proportional to exp {— ix 2 j (ujt) 2 }, 
in contrast to the typical case of random potential scattering when the amplitude of oscil- 
lations is proportional to exp {— it/(ujt)}. 

For the case of strong magnetic field (lot ^> 1) we obtain the DOS as a sequence of 
Gaussian broadened Landau levels with a width given by (|T2|). Note that the width does 
not increase with the magnetic field, like it does for the usual case of the potential scattering 
flT2|-|l4|. The reason for this unusual behavior is that the random phase acquired by the 
particle moving along its classical trajectory is proportional to the area enclosed by the 
trajectory rather than the length of the trajectory. 

To make this point clear consider the difference between our problem and that of the DOS 
oscillations in the uniform magnetic field in the presence of a smooth random potential V(r) 
with correlator (V(r)V(r')) = ^ 2 U exp{ — (r — r') 2 /^ 2 }, where £ ^> Ih is the correlation 
length. For our case the mean square random phase gained by the particle on its classical 
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trajectory with a winding number k is 

<0 2 > ~ (£f (h 2 )s no ~ (£f {h 2 )k 2 R 2 c ~ 

where t = 2irk/uj is the time spent by the particle on the trajectory. In case of potential 
scattering the analogous estimate gives instead: 

<0 2 > ~ ^R c k 2 ~ -^t 2 

where i>_p is the Fermi velocity. For both cases we can estimate r from the condition (0 2 ) ~ 1 
which leads in our case to the result (|T2|) and in the case of potential scattering reproduces 
results of [jl3| . 

Proportionality of (0 2 ) to the area, rather than to the length, of the trajectory actually 
determines the peculiarity of the random magnetic field scattering. Until the geometry 
of a certain problem is specified, any statement about single particle relaxation time will 
be senseless. We considered the DOS oscillations in a constant magnetic field, thus the 
characteristic area for this specific problem is the non-oriented area of the cyclotron circle: 
s n o = (vpi) 2 /An. If one would consider, say, the double-slit experiment QT5] the area which 



enters the problem will be determined by the characteristic length of the trajectories L and 
the characteristic angle a between them: s = L 2 sina. Thus one can estimate the single 
particle relaxation time in this case as follows: 1/r 2 ~ Emv^sina. Appearance of sin a in 
this expression reflects the dependence of 1/r on the geometry of the problem. 

Though the object proposed in [Kj (and later in ||) has, of course, the most evident 
choice of trajectory area for the two-point single particle Green's function, it is hard to 
imagine the physical problem in which it could enter. The Ref. addressed the question 
of the transport in the random magnetic field in terms of the path integral formalism. In 
this case, the relevant trajectories are of two-particle diffuson nature. Thus the the typical 
single particle trajectories, and their areas, considered in || are of dubious relevance to the 
transport problem. 

In conclusion, we discussed the single particle relaxation of 2D electrons in a spatially 
random magnetic field (in the presence of an additional uniform field). It was shown that 



any attempt to introduce a single particle relaxation time as a parameter of the Green's 
function of the type leads to an ambiguity since it violates the gauge invariance. The 
physically meaningful single particle relaxation time appears as a characteristic broadening 
of the Landau levels (jl2|) . Broadened Landau levels have the Gaussian shape fllTI). The 
width of the peaks does not increase with the magnetic field. The reason for this behavior 
is that the random phase acquired by the electron on a classical trajectory is proportional 
to the area enclosed by the trajectory, rather than to its length. 
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